A universal ionization threshold for strongly driven Rydberg states 
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We observe a universal ionization threshold for microwave driven one-electron Rydberg states of 
H, Li, Na, and Rb, in an ab initio numerical treatment without adjustable parameters. This sheds 
new light on old experimental data, and widens the scene for Anderson localization in light matter 
interaction. 
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The suppression of quantum transport across disor- 
dered media is one of the most spectacular consequences 
of destructive quantum interference. Originally predicted 
by Anderson [l| in his treatment of electrons propagating 
in disordered one dimensional lattices, Anderson localiza- 
tion has now become a general concept which prevails in 
abundant scenarios of coherent quantum transport in the 
presence of disorder 0, . Once it was realized that dy- 
namical chaos has the potential to substitute disorder in 
the long time evolution of low dimensional systems on the 
classical as well as on the quantum level :4] , it was natural 
to seek for a dynamical counterpart of Anderson localiza- 
tion in such a setting. Soon dynamical localization, i.e., 
the quantum suppression of diffusive energy growth in 
periodically driven, classically chaotic quantum systems 
was predicted 0, 0, f° r the periodically kicked rotor 
and for microwave driven hydrogen atoms initially pre- 
pared in highly excited Rydberg states. Whilst experi- 
ments on cold atoms Q - which closely mimic the kicked 
rotor in their center of mass motion - have confirmed 
this prediction in much detail, the situation remained 
controversial when it comes to the interpretation of ex- 
perimental results on Rydberg states of atomic h ydrogen 
and of alkaline atoms 

There, the presence of additional degrees of freedom, of 
the atomic continuum, and of a multielectron atomic core 
prevent a direct mapping onto the Anderson picture, and 
considerably complicate the unambiguous interpretation 
of experimental results. On the other hand, such atomic 
systems obviously extend the potential realm of dynam- 
ical localization - i.e. of Anderson localization on the 
energy axis - considerably, since they provide generic ex- 
amples of energy transport in light matter interaction. 

In this Letter we present a unifying theoretical picture 
which, for the first time, resolves all apparent inconsis- 
tencies between experimental observations on different 
atomic species under microwave driving 0, 0, 0, 
ITsl Tl9| . On the basis of ample numerical data obtained 
from an accurate, approximation-free treatment of the 
atomic excitation and ionization process, we identify a 
scaling law under which the ionization thresholds of dif- 
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ferent atomic species coalesce for sufficiently high driving 
field frequencies, independently of their element specific, 
unperturbed spectral structure. This universality of the 
ionization threshold provides strong support for the An- 
derson scenario as the underlying transport mechnism, 
clearly beyond all evidence provided so far. Furthermore, 
it allows a purely spectral interpretation of the observed 
localization mechanism, irrespective of the availability of 
a well defined classical analogue. 

Let us start with a condensed description of the phys- 
ical system under study: In a typical laboratory experi- 
ment 0, 0, 0, , atomic one electron Rydberg states 
prepared in an initial field free state po are exposed to 
a classical monochromatic microwave field of frequency 
uj and field strength F, for a duration t of typically sev- 
eral hundreds to several thousands of microwave periods 
T = 2ir/ui. After the atom field interaction, the per- 
centage P SU rv of atoms which did not ionize under the 
external driving is registered as a function of F and of 
the initial principal quantum number uq which labels po- 
At fixed uj and t, the ionization threshold field F(10%) 
which induces 1 — P surv = 10% is extracted from a se- 
ries of such measurements, over a large interval of no- 
For a microwave frequency which is of the order of the 
local energy spacing AE in the unperturbed Rydberg 
series, approx. N ~ 1/2uqAE photons have to be ab- 
sorbed by the atom to induce a transition from the ini- 
tial state po into the atomic continuum, i.e. to produce 
a nonvanishing ionization signal. For typical values of 
no ~ 25 . . . 100, this implies multiphoton transitions of 
the order 150 . . . 10, which can mediate efficient ioniza- 
tion only if they are composed of a sequence of near- 
resonant one photon transitions between Rydberg states 
separated by uj ± 8, with small detuning 8 -C uj. The 
increasing density of states within the Rydberg progres- 
sion towards the continuum threshold guarantees the ex- 
istence of such a sequence, provided uj is large enough 
to depopulate po via a first near resonant one photon 
transition. Under these premises, the intuitive analogy 
with the one dimensional Anderson model is established 
by identifying the sequence of near resonantly coupled 
Rydberg states with the neighbouring sites of a tight 
binding potential, with transition matrix elements essen- 
tially determined by the Rabi frequency which character- 
izes the one photon coupling. The one photon coupling 
constants are effectively randomized by the quasi ran- 
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dom distribution of detunings which is due to the non- 
linearity of the Rydberg pro gression to be set equal to 
m x w, with integer m |2(| . On the other hand, for 
AE ~ n^ 3 - which defines the spacing of the degener- 
ate energy levels of atomic hydrogen - the quantum me- 
chanical transport process along the energy axis can be 
compared to classically chaotic transport in phase space, 
since u> ~ n^ 3 defines strong nonlinear coupling between 
the unperturbed classical Kepler motion of the Rydberg 
electron and the external field, and hence leads to the 
efficient destruction of invariant tori in phase space for 
sufficiently large values of F 0, • Given this classical 
picture of strongly driven Coulomb dynamics, the results 
of typical experiments on atomic hydrogen are most of- 
ten represented in terms of scaled variables too — uj x Uq 
and Fo(10%) = F(10%) x Uq, what allows a direct iden- 
tification of the experimental (quantum) result with the 
associated classical phase space structure which is com- 
pletely determined by lu q and F Q . In such a plot, an 
increase of Fq(10%) is considered as a signature of dy- 
namical or Anderson localization, since classically diffu- 
sive/chaotic ionization leads to a monotonous decrease of 
Fn(10%) with ljq, for not too short atom-field interaction 
times [l8|. 

Despite the fact that no uniquely defined classical one 
particle dynamics is available for the excitation and ion- 
ization of nonhydrogenic initial states of alkaline Rydberg 
states, we compare in Fig. 1 the numerically generated 
scaled ionization thresholds of H, Li, Na, and Rb, as a 
function of ujq. The atoms are initially prepared in the 
low angular momentum state \hq£q = m a = 0), with 
£q the angular momentum and mo its projection on the 
polarization axis of the linearly polarized driving field. 
These results were obtained by diagonalization of the 
complex symmetric Floquet matrix which represents the 
complex dilated Floquet Hamiltonian of the driven atom, 
amended by the phase shift experienced by the electronic 
wave function of the Rydberg electron upon scattering 
from the multielectron core. Whilst details of our theo- 
retical and numerical treatment can be found elsewhere 
[23L |25| . let us just stress here that the atomic object 
under study is treated without any approximations on 
its dimensionality, that the atomic continuum is fully 
taken into account, and that no adjustable parameters 
are available. The results for the different species were 
obtained for fixed laboratory parameters w/27r = 36 GHz 
and t = 327 x T, and over a broad range 28 < no < 80 
of principal quantum numbers. Hence, representing the 
data in scaled units does not imply any a priori assump- 
tion on the relevant atomic energy scales, and the only 
possible cause of different ionization thresholds of differ- 
ent species is the element specific value of the nonvan- 
ishing quantum defects of the low angular momentum 
states. 

Clearly, we can identify three regimes [2^] in Fig. 1, 
which are defined by comparison of the hydrogen data 
with those for the nonhydrogenic alkaline initial states. 
In regime (I), at large scaled frequencies ljq > 1, all 



atomic species exhibit essentially identical ionization 
thresholds, irrespective of the specific structure of the un- 
perturbed Rydberg spectra characterized by quite vari- 
able quantum defects ranging from to 3.6, for £ < 3. 
Only small relative differences on top of the general trend 
of the ionization threshold in this frequency regime re- 
flect local spectral structures which differ from species 
to species and which determine the weight of individual 
atomic eigenstates in the field, in the decomposition of 
the initial field free state [ill 1231 [24J . These differences 
are particularly pronounced for the experimental hydro- 
gen data 0, |2| of the Stony Brook group also shown in 
this plot, since in these experiments (which, in addition, 
start from a statistical mixture po over the energy shell la- 
belled by no) the microwave is switched on rather slowly, 
allowing for a complicated sequence of (a)diabatic transi- 
tions during the switching period of the microwave field. 
Notwithstanding, the global agreement of numerical and 
experimental data in regime (I) is rather impressive. 

In regime (II), a strong discrepancy between the non- 
hydrogenic alkaline thresholds and those of atomic hydro- 
gen develops as we proceed to smaller scaled frequencies, 
i.e. to smaller pri ncip al quantum numbers at fixed labo- 
ratory frequency [l|| |21j . This sharp transition from (I) 
to (II) reflects the unavailability of a sequence of near 
resonant one photon transitions out of the hydrogen ini- 
tial state quite dramatically, since the hydrogen thresh- 
old can take values larger by more than almost one or- 
der of magnitude as compared to the alkaline threshold. 
Remarkably, the alkaline thresholds still coalesce in this 
interval of principal quantum numbers, and the increase 
of -Fb(10%) with increasing luq shows that the alkaline 
excitation and ionization process is still dominated by 
Anderson localization in this interval of scaled frequen- 
cies. 

Only in regime (III) is the energy of the driving photon 
too small to couple one photon transitions in the alka- 
line species, and the threshold saturates accordingly (as a 
matter of fact, it slowly increases with decreasing scaled 
frequency, what is however barely realized on the scale 
of the figure - see [2^] for an alternative presentation of 
the lithium data). 

In the inset in Fig. 1, we compare the above ioniza- 
tion threshold of Rb at ui/2tt = 36 GHz with the one 
at w/27r = 8.867 GHz, within the same loq interval and 
accordingly chosen values of no = 59 . . . 80, at the same 
interaction time t = 327 x T, and for a nonhydrogenic 
initial state with Iq = 1. Clearly, both thresholds conin- 
cide in scaled units, with a small systematic shift of the 
o;/27r = 8.867 GHz results towards lower values, what is 
simply due to the systematically larger values of no im- 
posed by the frequency scaling for the smaller laboratory 
value of the driving frequency [2^] ■ Obviously, the clas- 
sical scaling of frequency and amplitude cannot account 
for the semiclassical limit - which has to comply with 
the correspondence principle and hence must reproduce 
the classical ionization threshold which - by definition 
- is smaller than the quantum threshold in the Ander- 
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FIG. 1: Ionization thresholds of H (open diamonds, dotted line, and open circles, full line: laboratory experiments || l2l) : 
open circles, dash-dotted line: numerical experiments), Li (filled circles, full line), Na (open circles, dotted line), and Rb (filled 
diamonds, full line), in scaled units u>o and Fo(10%). All data were obtained for an interaction time t — 327 x T, at driving 
frequency lu/2tt = 36 GHz, with initial principal quantum numbers in the range no = 28 ... 80 except the ones represented 
by the open diamonds (here, lo/2tt — 9.95 GHz, and no = 32 . . . 90 (2^]). Hence, the comparison of alkali and hydrogen 
thresholds in scaled units does not imply any a priori assumption on the relvant atomic energy scales. Whilst all alkali data 
were obtained for nonhydrogenic initial states with angular momentum quantum numbers £q — and mo = 0, we observe a 
species-independent, universal ionization threshold in the high frequency range (I) - in strong support of the Anderson scenario 
as the dominant transport mechanism along the energy axis. Inset: Comparison of scaled Rb ionization thresholds, for driving 
frequencies u/2ir = 36 GHz, n = 38 . . . 66, £ = (filled triangles) and lu/2tt = 8.867 GHz, no = 59 . . . 80, £ = 1 (filled 
pyramids), respectively, at the same interaction time t = 327 x T. Clearly, the scaled thresholds essentially coincide, except 
for a systematic shift of the 8.867 GHz data to slightly smaller threshold values. This is a consequence of approaching the 
semiclassical limit h ~ rig — ► whilst keeping the scaled frequency fixed |2ll |22|. Also shown are experimental data (open 
circles) obtained at 8.867 GHz, though for much longer interaction times t = 44335 xT 13]. Whilst the general uio dependence 
is identical (notably the transition between frequency ran ges (HI) and (II)) for numerical and experimental data, the typically 
algebraic time dependence of the ionization threshold (l3l |23|| induces smaller thresholds for longer interaction times fl9ll24|. 



son/dynamically localized regime. 

The inset also contains experimental ionization thresh- 
olds of Rb £o — 1 states under lu/2tt = 8.867 GHz driving, 
though for much longer interaction times t — 44335 x T 
[l3[ . Clearly, due to the nonvanishing decay rate of the 
eigenstates of the atom in the field, the experimentally 
observed ionization threshold is smaller than the numer- 
ical one obtained for shorter interaction times, but the 
transition between regimes (II) and (III) occurs in the 
same range at the same of loq [23|! Furthermore, we see 
that the time dependence of the threshold is stronger for 
larger frequencies, thus slightly affecting the slope of the 
threshold curve - an observation also made in labora- 



tory experiments on Li [T^. Indeed, a detailed study of 
the time dependence of the ionization threshold reveals 
a generic, algebraic decay law F (10%) ~ t~ 7 , where 7 
tends to increase slowly with uj$ [2jj. Whilst we could 
not reliably access the experimental regime of interaction 
times, since this requires extremely large a basis set such 
as to achieve numerical convergence of the decay rates 
better than 10 -13 a.u., typical values of 7 are consistent 
with the gap opening up between the numerical and ex- 
perimental 8.867 GHz data in Fig. 1 (note, however, that 
the precise value of 7, at a given value of uq and u>, will 
depend on the local spectral structure and on the partic- 
ular envelope of the microwave pulse experienced by the 
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atoms |l]J, |24J , which our numerical approach does not 
account for). 

Given our observations in Fig. 1, the following picture 
emerges: Atomic one electron Rydberg states exhibit a 
universal ionization threshold if the tight binding pic- 
ture sketched initially can be realized by a sequence of 
near resonant one photon transitions connecting the ini- 
tial field free state to the atomic continuum [20]. Con- 
sequently, in regime (I), for u>q > 1, i.e. for photon en- 
ergies larger than the hydrogcnic level spacing, alkaline 
atoms and atomic hydrogen exhibit identical thresholds 
and ignore the additional nonhydrogenic level structure 
provided by the nonvanishing quantum defects of the dif- 
ferent species. Also for alkaline atoms the scaling law 
inherited from the classical Coulomb dynamics prevails, 
since the energy spacing between the nonhydrogenic ini- 
tial state and the hydrogenic manifold exhibits the same 
dependence on the principal quantum number as the hy- 
drogenic energy splitting |2^. Nonhydrogenic states of 
alkali atoms exhibit Anderson/dynamical localization in 
regime (II), in the absence of well defined classical one 
electron dynamics, simply due to the availability of a se- 
quence of near resonant one photon transitions to the 
continuum. Note that in this regime this spectral sub- 



structure is provided by the scattering of the electron 
from the atomic multielectron core which, on the scale of 
a typical Rydberg orbit in the energy range considered 
here, acts as a point scatterer [2(|. Given the tremendous 
enhancement of the alkaline ionization yield as compared 
to atomic hydrogen in this parameter range, this rep- 
resents the arguably most dramatic, directly observable 
signature of core scattering so far observed in a chaotic 
quantum system: here, core scattering restablishes An- 
derson localization! As compared to the core scattering 
effect the nonvanishing time dependence of the ionization 
threshold is but a correction: even increasing the interac- 
tion time by a factor 10 . . . 100 does lower the threshold 
only weakly. All published experimental data which were 
so far interpreted in terms of dynamical localization were 
clearly obtained in regimes (II) and (III). The universal 
threshold in regime (I) was a theoretical prediction [23j . 
which, however, has been confirmed in very recent labora- 
tory experiments [III on one electron Rydberg states 
of Li and Sr. 

Support as a Grand Challenge project at the Leibniz- 
Rechenzentrum of the Bavarian Academy of Sciences is 
most gratefully acknowledged. 



[1] P. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] D. S. Wiersma, P. Bartolini, A. Lagendijk, and R. Righ- 

ini, Nature (London) 390, 671 (1997). 
[3] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 

(1993) . 

[4] O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. 

Lett. 52, 1 (1984). 
[5] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. 

Rev. Lett. 49, 509 (1982). 
[6] G. Casati, B. V. Chirikov, and D. L. Shepelyansky, Phys. 

Rev. Lett. 53, 2525 (1984). 
[7] R. Bliimel and U. Smilansky, Phys. Rev. Lett. 52, 137 

(1984). 

[8] F. L. Moore, J. C. Robinson, C. Bharucha, P. E. 
Williams, and M. G. Raizen, Phys. Rev. Lett. 73, 2974 

(1994) . 

[9] E. J. Galvez, B. E. Sauer, L. Moorman, P. M. Koch, and 
D. Richards, Phys. Rev. Lett. 61, 2011 (1988). 

[10] J. E. Bayfield, G. Casati, I. Guarneri, and D. W. Sokol, 
Phys. Rev. Lett. 63, 364 (1989). 

[11] H. P. Breuer, K. Dietz, and M. Holthaus, J. Phys. B 22, 
3187 (1989). 

[12] P. Fu, T. J. Scholz, J. M. Hettema, and T. F. Gallagher, 

Phys. Rev. Lett. 64, 511 (1990). 
[13] M. Arndt, A. Buchleitner, R. N. Mantegna, and 

H. Walther, Phys. Rev. Lett. 67, 2435 (1991). 
[14] R. Graham, Comm. At. Mol. Phys. 25, 219 (1991). 
[15] T. F. Gallagher, Comm. At. Mol. Phys. 25, 159 (1991). 
[16] P. Koch, L. Moorman, and B. Sauer, Comm. At. Mol. 

Phys. 25, 165 (1991). 
[17] S. Fishman and D. L. Shepelyansky, Europhys. Lett. 16, 

643 (1991). 

[18] O. Benson, A. Buchleitner, G. Raithel, M. Arndt, R. N. 
Mantegna, and H. Walther, Phys. Rev. A51, 4862 



(1995) . 

[19] M. W. Noel, W. M. Griffith, and T. F. Gallagher, Phys. 

Rev. A 62, 63401 (2000). 
[20] N. Brenner and S. Fishman, Phys. Rev. Lett. 77, 3763 

(1996) . 

[21] P. M. Koch and K. A. H. van Leeuwen, Phys. Rep. 255, 
289 (1995). 

[22] A. Buchleitner, D. Delande, and J. C. Gay, J. Opt. Soc. 
Am. B 12, 505 (1995). 

[23] A. Krug, Ph.D. thesis, Ludwig-Maximilians- 
Universitat Miinchen (2001), http://edoc.ub.uni- 
muenchen.de / archive/00000336 / . 

[24] A. Buchleitner, D. Delande, J. Zakrzewski, R. N. Man- 
tegna, M. Arndt, and H. Walther, Phys. Rev. Lett. 75, 
3818 (1995). 

[25] A. Krug and A. Buchleitner, Phys. Rev. A 66, 53416 
(2002). 

[26] T. Jonckheere, B. Gremaud, and D. Delande, Phys. Rev. 
Lett. 81, 2442 (1998). 

[27] T. Gallagher, oral contribution, 296. WE-Heraeus- 
Seminar "Chaos and Quantum Transport" , Bad Honnef, 
23-27 March 2003. 

[28] H. Maeda and T. Gallagher (2004), preprint. 

[29] As a matter of fact, this scaling law also allows to un- 
derstand experimental data on Rb at ui/2n — 12.6 GHz 
[Tij|. where a prominent enhancement of the ionization 
threshold was observed at no = 89, though could not 
be explained: In Coulomb scaling, this corresponds to 
a scaled frequency luo — 1-3, and the local maximum of 
fb(10%) thus is a remnant of the principal nonlinear res- 
onance island in the classical phase space of the driven 
two-body system [2l|. 



